Spin-state transition and phase separation in multi-orbital Hubbard model 
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We study spin-state transition and pliase separation involving tliis transition based on the milti-orbital Hubbard 
model. Multiple spin states are realized by changing the energy separation between the two orbitals and the on- 
site Hund coupling. By utilizing the variational Monte-Carlo simulation, we analyze the electronic and magnetic 
structures in hole doped and undoped states. Electronic phase separation occurs between the low-spin band 
insulating state and the high-spin ferromagnetic metallic one. Difference of the band widths in the two orbitals 
is of prime importance for the spin-state transition and the phase separation. 
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Novel electric and magnetic phenomena observed in cor- 
related electron systems are responsible for competition and 
cooperation between multi-electronic phases with delicate en- 
ergy balance. These are owing to the internal degrees of free- 
dom of electrons, i.e. spin, charge and orbital, under strong 
electron correlation, and their coupling with crystal lattice.-iS 
In some transition-metal ions, there is an additional degree of 
freedom, termed the spin-state degree of freedom, i.e. mul- 
tiple spin states due to the different electron configurations 
in a single ion. One prototypical example is the perovskite 
cobaltites Ri^^^AxCoOi {R: rare earth ion. A: alkaline earth 
ion) where transitions between the multiple spin states occur 
by changing carrier concentration, temperature and so on. In 
Co^+ with the configuration, there are three possible spin 
states, the high-spin (HS) state (£^12^) with an amplitude of 
5 = 2, the intermediate-spin (IS) one (e^fl^) with 5=1, and 
the low-spin (LS) one (tig) with 5 = 0. 

Several magnetic, electric and transport measurements have 
been carried out in the insulating and metallic cobaltites. It 
is known that LaCoOs is a non-magnetic LS band-insulator 
(BI) at low temperatures, although there is still controversy 
in the spin-state transition and the IS state at finite temper- 
ature .^i^'^'^'^ In high hole doping region of x > 0.3 — 0.4 in 
Lai-^SrvCoOs, the ferromagnetic (FM) metallic state was ex- 
perimentally confirmed. In the lightly hole doped region be- 
tween the two, a number of inhomogeneous features in mag- 
netic, electric and lattice structures have been reported ex- 
perimentally. Spatial segregation of hole-rich FM regions 
and hole-poor insulating ones have been suggested by the 
neutron diffraction, the electron microscopy, NMR and so 
oni2i2ii2iii Magnetic/non-magnetic clusters have been found 
by the small-angle and inelastic neutron scattering experi- 
mentsji^ii^iii It is widely believed that the observed giant 
magneto-resistance effect in the lightly doped region results 
from the electronic and magnetic inhomogeneityJ^ 

Electronic phase separation (PS) phenomena in transition- 
metal compounds have been studied extensively and inten- 
sively, in particular, in the high Tc superconducting cuprates 
and the colossal magnetoresistive manganitesi '^i^^-^^''^ In 
these materials, the long-range spin/orbital orders in the Mott 
insulating phases and their melting by carrier doping are of 
essence in the electronic PS. The exchange energy for the lo- 



calized spins/orbitals and the kinetic one for the itinerant elec- 
trons are gained in spatially separate regions. On the other 
hand, in the present case, the non-magnetic band insulator is 
realized in the insulating phase, and the spin-state transition is 
brought about by carrier doping. Thus, the present phenom- 
ena belong to a new class of the electronic PS in correlated 
system, although only a little theoretical studies have been 
done until now. In this paper, we address the issues of the 
spin-state transition and the PS associated with this transition 
by analyzing the multi-orbital Hubbard model. We examine 
the electronic structures in hole doped and undoped systems 
by utilizing the variational Monte-Carlo (VMC) method. We 
find that, between the non-magnetic BI and the HS FM metal, 
the electronic PS is realized. We claim that the different band 
widths play an essential role in the present electronic PS. 

We set up a minimal model, the two-orbital Hubbard 
model j'VV^^^ where the spin-state degrees of freedom and 
a transition between them are able to be examined. In each 
site in a crystal lattice, we introduce two orbitals, termed A 
and B, which represent one of the eg and t2g orbitals, respec- 
tively. Anisotropic shape of the orbital wave function is not 
concerned. An energy difference between the two orbitals is 
denoted by A = — eg > where Ca (Eb) is the level energy 
for A(B). When the electron number per site is two, the lowest 
two electronic states in a single site are \B^) and IA'b') with 
triplet spin state which are termed the LS and HS states in the 
present model, respectively. The explicit form of the model 
Hamiltonian is given by 

ia {ij)ya 

iy iaa' 

^ H ^lAa^iBac]g^,CiAa' - fJ^cJ^Ciy^C^^^Ciyi, (1) 
iaa' iy 

where Ciya is the annihilation operator of an electron at site / 
with orbital /(= A,B) and spin CJ(=t, i), and riiya = c^fjCiya 
is the number operator A subscript f takes A{B), when 7 is 
B{A). We assume that the transfer integral is diagonal with re- 
spect to the orbitals and jf^l > Ifg], both of which are justified 
in perovskite cobaltites. In most of the numerical calculations, 
a relation tg/tA = 1/4 is chosen. As the intra-site electron in- 
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teractions, we introduce the intra- and inter-orbital Coulomb 
interactions, U and U' , respectively, the Hund coupling / and 
the pair-hopping J' . The relations U = U' + 2J and J ^J' sat- 
isfied in an isolated ion are assumed. In addition, we introduce 
the relation U = 4/ in the numerical calculation. 

We adopt the VMC method where the electron correlation 
is treated in an unbiased manner and simulations in a large 
cluster size are possible. For simplicity and a limitation in the 
computer resource, we introduce two-dimensional square lat- 
tices with a system size of N = Lr (L < 6) and the periodic 
and anti-periodic boundary conditions. The number of elec- 
tron is Ne, and the hole concentration per site measured from 
Ne = 2N is denoted as x = {2N — Ne)/N. The variational wave 
function is given as a product form of = G|<I>) where G is 
the correlation factor and |<I>) is the one-body wave function. 
The two types of the wave function are considered in |<t>): 
the Slater determinant obtained by the second term in Eq. ([T]), 
and that for the HS antiferromagnetic (AFM) order given by 
applying the Hartree-Fock approximation to the third term in 
Eq. ([Til. In the latter, the AFM order parameter is treated as a 
variational parameter. We assume the Gutzwiller-type correla- 
tion factor n,7(l — ^i^ii) where / indicates the local electron 
configurations, is the projection operator at site / for the 
configuration I, and £,i is the variational parameter. Here we 
introduce the 10 variational parameters for the 10 inequivalent 
electron configurations in a single site^— The fixed-sampling 
method is used to optimize the variational parametersi^l In 
addition to the standard VMC method, we improve the vari- 
ational wave function by estimating analytically the weights 
for the configurations which are sampled by the MC simu- 
lations. This method is valid for the LS state and reduces the 
CPU time by more than one order. In most of the calculations, 
10"*^ — 10^ MC samples are adopted for measurements. 

We start from the case at x = where the average 
electron number per site is two. The electronic states 
obtained by the simulation are monitored by the total 
spin amphtude defined by = il/N)Y,i{Sf) where S,- = 
Y.ySiy = (l/2)^jj/ycjy^(7.„/c,yy is the spin operator with the 
Pauli matrices a, the spin correlation function Sy{q) = 
(4/A^)^,-,e"i ('"'"'"j^(5-y5j^), and the momentum-distribution 
function «y(k) = (l/2)^^(c[^^Cky(7) where Ckya is the 

Fourier transform of Cjya- Size dependences of and ^^(q) 
in L = 4 — 8 are within a few percent. We obtain the three 
phases, the HS Mott insulator (MI), the LS BI and the metal- 
lic (ML) phase. In the HS-MI phase, is about L6 being 
about 80% of the maximum value for S = I. A sharp peak in 
^^(q) at q = {7t,n) and no discontinuity in ny(k) imply that 
this is the AFM MI. In the LS-BI phase, «A(k) [«B(k)] is al- 
most zero (one) in all momenta, and ~ 0. In the ML phase, 
discontinuous jumps are observed in both «/i(k) and «B(k). 
The electron (hole) fermi surface is located around k — (0,0) 
[{n, n)] in the A (B) band; this is a semi metal. A value of 
is about 0.3, and no remarkable structure is seen in ^^(q). 

The phase diagram at x = is presented in Fig.[T] The error 
bars imply the upper and lower bounds of the phase bound- 
ary, and symbols are plotted at the middle of the bars. In the 
region of large A (7), the LS-BI (HS-MI) phase is realized. 
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FIG. 1 : (color online) Phase diagrams at x = 0. A ratio of the electron 
transfers is taken to be = 1/4 in (a) and ts/tA = I in (b). In 
(b), filled squares and open circles are for the results obtained by the 
VMC method and the previous DMFT one in Ref.[j3, respectively. 
Broken curves are guides for eyes. Stars represent the parameters 
where the carrier dopings are examined. 



and between the two with small A and /, the ML phase ap- 
pears. To compare the present results with the previous ones 
calculated by the dynamical-mean field theory (DMFT)fi^ we 
present, in Fig.[TJb), the phase diagram where the two transfer 
integrals are chosen to be equal, i.e. Ib/Ia = 1. Although the 
global features in the phase diagrams are the same with each 
other, the HS-MI phase obtained by the VMC method appears 
in a broader parameter region than that in DMFT, in partic- 
ular, near the boundary of the HS-MI and ML phases. This 
is because the AFM long-range order in the HS-MI phase is 
treated properly in the VMC method. We have confirmed that 
the phase boundaries obtained by the VMC method where the 
AFM order is not considered almost reproduce the DMFT re- 
sults. 

Now we show the results at finite x. Holes are introduced 
into the LS-BI phase near the phase boundary with the pa- 
rameter values of (A/fA,//fA) = (12.2,4) and (8.25,2.5) [see 
Fig.[T]. By changing the initial conditions in the VMC simu- 
lation, we obtain the following four states: i) the LS-ML state 
where «A(k) is almost zero in all k, and the fermi surface is 
located in the B band around k = {n,7z), ii) the FM HS-ML 
state where nB(k) is about 1/2 in all k, the fermi surface is 
in the A band, and Sy{q) has a sharp peak at q = (0,0), iii) 
the AFM HS-ML state where the fermi surface exists in the A 
band around k ~ (?r,0), and ^^(q) has a peak at q = {Tt,n), 
and iv) the mixed state where the wave function is a linear- 
combination of the LS-ML and FM HS-ML states. 

In Fig.|2ja), the energy expectation values E = (J^) for the 
several states in {A/tA,J /Ia) ~ (12.2,4) are plotted as func- 
tions of X. The transfer integrals are chosen to be /Ia = 1/4. 
To show the numerical data clearly, we plot E' ^ {E /Ia) + Cx 
with a numerical constant C, instead of E. This transforma- 
tion does not affect the Maxwell's construction introduced be- 
low. The results in the AFM HS-ML are not plotted, because 
of their higher energy values than others. We also present, 
in Fig. [3] a ratio of the LS sites to the LS and HS ones in 
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FIG. 2: (color online) Hole concentration dependences of the energy 
expectations for several states at {A/tA,J/tA) = (12.2,4) in (a), and 
those at {A/tf^,J /tp^) = (8.25,2.5) in (b). Broken lines are given by 
the Maxwell's construction. A ratio of the electron transfers is taken 
to be tB jtA = 1/4. A constant parameter C in the definition of E' is 
taken to be 8.2 in (a) and 5.25 in (b). 



the mixed states defined by = nis/ [uis + nHs)- Here tils 
(tiHs) is a number of the sites where the LS (HS) state is re- 
alized. As shown in Fig. |2a), the LS state, where holes are 
doped into the B band, is destabilized monotonically with in- 
creasing X. On the other side, in a region of x > 0.5, the 
FM HS-ML state is realized. In between the two regions, 
the mixed state is the lowest energy state. The mixed state is 
smoothly connected to the LS and HS ones in the low and high 
X regions, respectively. As shown in Fig. [3] a discontinuous 
jump in the mixed state is seen around x ~ 0.25; the system is 
changed from the LS dominant mixed state into the HS dom- 
inant one with x. It is noticeable that the E' versus x curve in 
the mixed state is convex in the region of < x < 0.33. That is, 
by following the Maxwell's construction, the PS of the LS-BI 
and the FM HS dominant mixed states is more stabilized than 
the homogeneous phase in this region of x. In the Fig. |2jb), 
we show the results in (A/fA,//fA) ~ (8.25,2.5) where the 
system at jc = is close to the ML phase [see Fig.lTJa)]. The 
PS appears, but its region is shrunken. 

The magnetization per site in the lowest energy state de- 
fined by M{x) = (1/A^)(I;5;') is plotted in Fig. [3 A zero 
magnetization at x = reflects the LS-BI ground state. In a 
high doped region of x > 0.33, the magnetization data almost 
follow a relation M(x) ~ (1 +x)/2; the system is expected 
to consist of the A^/2 HS sites, the (1/2 —x)N LS ones, and 




FIG. 3: (color online) A ratio of the LS sites to the LS and HS 
ones in the mixed state, Ris, and magnetization M(x) as func- 
tions of the hole concentration x. A broken line connecting data at 
M(x = 0) and M{x = 0.33) is drawn by the Maxwell's rule. For 
comparison, we plot a M{x) = x/2 curve which is expected from 
the hole doping in the LS-BI phase. Parameters are chosen to be 
{A/lA,J/tA) = (12.2,4) and/B/fA = 1/4. 




FIG. 4: (color online) Hole concentration dependences of the energy 
expectations for several states where the electron transfer integrals 
are chosen to be equal as ts /t/^ = I. Other parameters are taken to be 
(A/tA,J /tA) = (12.2,4), and a constant parameter C in the definition 
of E' is taken to be 8. 



the xA^ singly electron occupied ones. In this scheme, we ob- 
tain Rls = ( 1 ~ 2x) / (2 — 2x) which is consistent with the nu- 
merical data of Rls in x > 0.33. Between x = and 0.33, 
where the PS is realized, M(x = 0) and M(x = 0.33) are con- 
nected by a straight line according to the volume-fraction rule 
in the Maxwell's construction. The slope of M(x) is about 
three times higher than M(x) = x/2 which is expected in the 
hole doping into the LS-BI phase. This is qualitatively consis- 
tent with the experimental observations in the magnetization 
where doped holes induce high spin valucj^^ 

We now address an origin of the electronic PS where the 
spin-state degree of freedom is concerned. In Fig. |4] we 
present the hole concentration dependence of the energy ex- 
pectations where the band widths are set to be equal with each 
other, ffi/f/i = 1 . As well as the calculation in Fig. |2a), the en- 
ergy parameters are taken to be {A/tA,J /(a) = (12.2,4) which 
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FIG. 5: (color online) Schematic density-of-states in the LS-BI state 
at X = and that in the HS-ML one in a high hole doped region. 



is close to the LS-HS phase boundary at jc = [see Fig.lTJb)]. 
The mixed state is not obtained in the simulation. In all region 
of X up to X = 0.45, the LS state is the lowest ground state, and 
neither the spin-state transition nor the PS occur. The differ- 
ence of the band widths in the two orbitals is of essence in the 
electronic PS phenomena. 

To clarify the mechanism of PS furthermore, schematic pic- 
tures of the density of states (DOS) in the LS-BI at x = and 
the FM HS-ML in a high hole doped region are presented in 
Fig. ID For simplicity, detailed shapes of DOS are not taken 
into account. In LS-BI state at x = 0, the fermi level is lo- 
cated inside of the band gap between the A and B bands. The 
band width in the A band is larger than that in B. On the other 
hand, in the FM HS-ML state which is realized in x > 0.5 in 
Fig. 121 a), the system is a doped MI with ferromagnetic spin 
polarization. The fermi level is located in the A band. Be- 
cause of the large band width in the A band, there is a large 
kinetic energy gain in comparison with the doped LS-BI state 
where the fermi level is located in the B band in the rigid band 
scheme. This kinetic energy gain is the origin of the spin state 
transition by doping. It is shown in Fig.|4]that, when the equal 
band widths are assumed, the E' v.s. x curves for the LS-ML 
and FM HS-ML states are almost parallel and do not cross 



with each other. This data implies that there is no difference 
in the kinetic energy gains for the two states, when the band 
widths are assumed to be equal. The present PS phenomena 
are also attributed to this band width difference as follows. In 
the rigid-band sense, by doping of holes in the LS-BI state, 
the fermi level falls into the top of the B band from the middle 
of the gap in Fig. |3a). If we suppose that this state is real- 
ized in a low x region and is transferred into the FM HS-ML 
state shown in Fig. |5jb) with increasing x, the fermi level is 
increased with increasing hole concentration because of the 
different band widths. This is nothing but the negative charge 
compressibility K ~ {djJ./ dx) < with the chemical potential 
jj., i.e. appearance of the electronic PS. 

Finally, we discuss implications of the perovskite 
cobaltites. The obtained PS between the insulating nonmag- 
netic state and the hole-rich FM one is qualitatively consis- 
tent with the inhomogeneity suggested by a number of ex- 
periments. The PS and the spin-state transition are attributed 
to the band-width difference of the two bands correspond- 
ing to the Cg and t2g bands in the perovskite cobaltites. This 
electronic PS is robust by changing the model parameter val- 
ues, except for ts/tA, when the non-doped system is located 
near the phase boundary between the LS-BI and HS-ML The 
present phenomena are different from the previous PS's dis- 
cussed in the high-Tc cuprates and the manganites where the 
long-range spin/orbital orders are realized in the Mi's; the 
spatial segregations occur between the long-range ordered MI 
and the ML states where the superexchange interaction en- 
ergy and the kinetic one of doped holes are separately gained 
in the different spatial regions. Our scenario of the PS based 
on the band-width difference may be checked experimentally 
by adjusting the tolerance factor, i.e. the Co-O-Co bond an- 
gle; the smaller tolerance factor implies the smaller (larger) 
band width in the eg {t2g) orbitals, and suppression of the PS. 
Detailed values of x where the PS is realized, and a typical 
size of the clusters remain as questions. Several factors not 
considered here, the intermediate-spin state, the long-range 
Coulomb interaction, the lattice volume depending on the spin 
states, and so on, are required to answer these questions. 

Authors would like to thank H. Yokoyama and 
H. Takashima for their valuable discussions. This work 
was supported by JSPS KAKENHI, TOKUTEI from 
MEXT, and Grand challenges in next-generation integrated 
nanoscience. 



S. Maekawa, et al. Physics of Transition Metal Oxides, (Springer 
Verlag, Berlin, 2004), and references therein. 

2 M. Imada, et al. Rev. Mod. Phys. 70, 1039 (1998). 

3 S. Yamaguchi et al. Phys. Rev. B 53, R2926 (1996). 

4 M. A. Korotin, et al. Phys. Rev. B 54, 5309 (1996). 

5 M. W. Haverkort, et al. Phys. Rev. Lett. 97, 176405 (2006). 
^ S. Noguchi, et al. Phys. Rev. B 66, 094404 (2002). 

^ Y. Kobayashi, et al. Phys. Rev. B 72, 174405 (2005). 
^ M. Itoh, et al. J. Phys. Soc. Jpn. 63, I486 (1994). 

R. Caciuffo, etal. Phys. Rev. B 59, 1068 (1999). 
'° P L. Kuhns, etal. Phys. Rev. Lett. 91, 127202 (2003). 



■ ' A. Ghoshray, et al. Phys. Rev. B 69, 064424 (2004). 
'2 J. Wu, etal. Phys. Rev. Lett. 94, 037201 (2005). 
13 D. Phelan, etal. Phys. Rev. Lett. 96, 027201 (2006). 
'4 D. Phelan, etal. Phys. Rev. Lett. 97, 235501 (2006). 
■5 E. L. Nagaev, Phys. State. Sol. (b) 186, 9 (1994). 

E. Dagotto, The Physics of Manganites and Related Compounds, 

(Springer- Verlag, Berlin 2003). 
1^ S. Okamoto, etal. Phys. Rev. B 61, 451 (2000). 

K. I. Kugel, etal. Phys. Rev. Lett. 95, 267210 (2005) 

P Werner, and A. J. Millis, Phys. Rev. Lett. 99, 126405 (2007). 
2" K. Sano and Y. Ono, J. Phys. Soc. Jpn. 72, 1847 (2003). 



5 



K. Kobayashi, and H. Yokoyama, Physica C 445-448, 162 (2006). 

K. Kubo, Phys. Rev. B 79, 020407 (2009). 

The 10 inequivalent configurations considered here are |0), \A'^), 

\B"), lA^B^), lA'^S"), IflTgi), \aU^), \A"B^B^), 1aU1s">, 

lAUJ-fiTgi) with <T =t and i, where A'^{B'^) impies that the A 

(B) orbital is occupied by the up (down) spin electron. 

C. J. Umrigar, et al. Phys. Rev. Lett 60, 1719 (1988). 



25 J. Okamoto, et al. Phys. Rev. B 62, 4455 (2000). 

^ Present address; The bank of Tokyo-Mitsubishi UFJ, Tokyo, 
Japan. 

* Present address: Chiba Institute of Technology, Tsudanuma, 
Chiba 275-0016, Japan. 



